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We study the path integral formulation of Friedmann universe iilled with a massless scalar field in loop quan- 
tum cosmology. All the isotropic models of A: = 0, +1, —1 are considered. Since the transition amplitude in the 
deparameterized framework can be expressed in terms of group averaging, the path integrals can be formulated 
for both deparameterized and timeless frameworks. It turns out that the effective Hamiltonian derived from 
the path integral in deparameterized framework is equivalent to the effective Hamiltonian constraint derived 
from the path integral in timeless framework, since they lead to same equations of motion. Moreover, the ef- 
fective Hamiltonian constraints of above models derived in canonical theory are confirmed by the path integral 
formulation. 
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^ . I. INTRODUCTION 



A core concept of general relativity is that physical theory should be background independent, which indicates that the gravity 
and matter evolve as a whole entity without anything pre-existed. This idea is being carried out in loop quantum gravity 
(LQG)[l-4], which is a non-perturbative and background independent quantization of general relativity. The construction of 
^ ' LQG inspired the research on spin foam models (SFMs) [5], which are proposed as the path-integral formulation for LQG. 
y^. In SFMs, the transition amplitude between physical quantum states is formulated as a sum over histories of all physically 
5—1 ' appropriate states. The heuristic picture is the following. One considers a 4-dimensional spacetime region M bounded by two 
spacelike 3-surfaces, which would correspond to the initial and final states for the path integral. The feature of SFMs is to 
employ the spin network states in LQG to represent the initial and final states. A certain path between the two boundaries is a 
Tj" quantum 4-geometry interpolated between the two spin networks. The interpolated spin foams can be constructed by considering 
the dual triangulation of M and coloring its surfaces with half integers and edges with suitable intertwiners. In order to obtain 
the physical inner product between the two boundary states, we have to sum over all possible triangulations and colorings [1]. 
This is certainly difficult to be achieved. 

On the other hand, the idea and technique of LQG is successfully carried on in the symmetry- reduced models, known as loop 
quantum cosmology (LQC) [6]. Since the infinite degrees of freedom of gravity have been reduced to finite ones by homogeneity 
(and isotropy), LQC provides a simple arena to test the ideas and constructions of the full LQG. Therefore, it is also desirable 
to test the idea and construction of SFMs by LQC models [7]. It was shown in Refs.[8, 9] that the transition amplitude in the 
■ deparameterized LQC of fc = Friedmann universe equals the physical inner product in the timeless framework, and concrete 
• • [ evidence has been provided in support of the general paradigm underlying SFMs through the lens of LQC. Recently, the effective 
. !^ ' Hamiltonian constraint of fc = LQC was also derived from the path integral formulation in timeless framework in Ref. [10], and 
[ the effective measure introduces some conceptual subtleties in arriving at the WKB approximation. The purpose of this paper 
;h I is to study the path integral formulation of LQC mimicking the construction of SFMs and following [7-10]. Besides the fc = 
. Priedmann-Rorberston- Walker (FRW) model, the fc = ±1 models will also be considered. It turns out that the physical inner 

product between boundary states, defined by group averaging approach, can be interpreted as sum over all possible quantum 
states. Moreover, we will be able to derive the effective action and thus Hamiltonian of the above LQC models by the path 
integral formalism in deparameterized framework as well as timeless framework. The effective Hamiltonian coincide with those 
obtained by canonical approach. Hence the path integral and canonical formulations of LQC confirm each other The effective 
Hamiltonian constraints of above models in canonical theory are also confirmed by the path integral formulation. Starting from 
the canonical LQC, we will derive the path integral formulation and the effective Hamiltonian of fc = 0, +1, —1 FRW models 
respectively in the following Sections 2, 3, and 4. To simplify the calculations, we will employ the simplified LQC treatments 
[11]. 
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II. PATH INTEGRAL FORMULATION OF fc = LQC 

The model that we are considering is the spatially flat (fc = 0) FRW universe filled with a massless scalar field. In the 
kinematical setting, one has to introduce an elementary cell V in the spatial manifold and restricts all integrations to this cell. 
Fix a fiducial flat metric "g^^ and denote by Vo the volume of V in this geometry. The gravitational phase space variables — the 

connections and the density weighted triads — can be expressed as ^ cVo ^ ' ' °< and Ef ^ pVo ^ ' ' °ef, where 
, °ef) are a set of orthonormal co-triads and triads compatible with °q^f^ and adapted to V. p is related to the scale factor a 

via \p\ = Vo^^a^- The fundamental Poisson bracket is given by: {c, p} ~ SnG^/i, where G is the Newton's constant and 7 
the Barbero-Immirzi parameter The gravitational part of the Hamiltonian constraint reads Cgrav = ~6c^y^|p[/7^. To apply the 
area gap A of full LQG to LQC [12], it is convenient to introduce variable fi satisfying 

fl^\p\ = A=(4VTi7rj)£l, (2.1) 

where = Gh, and new conjugate variables [13, 14]: 

^;:=2^/3^^-^ b -.^ flc. (2.2) 

The new canonical pair satisfies {6, v} = j-. On the other hand, the matter phase space consists of canonical variables and 
P0 which satisfy {t/jjp^} — 1. To mimic LQG, the polymer-like representation is employed to quantize the gravity sector The 
kinematical Hilbert space for gravity then reads "H^™" = i^(RBohi-, c^/^h), where RBohr is the Bhor compactification of the 
real line and diiH is the Haar measure on it [15]. It turns out that the eigenstates of volume operator v, which are labeled by 
real number v, constitute an orthonomal basis in Hf,™" ^s {vi\v2) = Svi.v2- For the scalar matter sector, one just uses the 
standard Schrodinger representation for its quantization, where the kinematical Hilbert space is W^^^^^ = i^(]R, d(j)). The total 
kinematical Hilbert space of the system is a tensor product, T-Lkm = '^ki'n ® ^Sn**' o^ above two Hilbert spaces. As a 
totally constrained system, the dynamics of this model is reflected in the Hamiltonian constraint Cgrav + Cmatt = 0. Quantum 
mechanically, physical states are those satisfying quantum constraint equation 

(Cgrav + a,att)*(«, =0. (2.3) 

It is not difficult to write the above equation as [16] 



C-^{v,cj>) = {?^-%)^{v,c^)^0, (2.4) 



where 8 is a self-adjoint and positive-definite operator On one hand, Eq.(2.4) indicates that we can get physical states by group 
averaging kinematical states as [9] 

/oo ^ 
da e'"^2|p^| (2.5) 
-OO 

and thus the physical inner product of two states reads 

/•oo 

( f\9 )pky = {^f\g )= da {fW''^2\p^\ \g). (2.6) 



On the other hand, Eq.(2.4) can also be written as 

dl^{v,(j)) + m{v,(t))=Q. (2.7) 

The similarity between Eq.(2.7) and Klein-Gorden equation suggests that one can regard as internal time, with respect to which 
gravitational field evolves. In this deparameterized framework, we focus on positive frequency solutions, i.e., those satisfying 

(j,) = Ve*+(w, = (j)). (2.8) 

The transition amplitude in deparameterized framework is then given by 

Aivf.cj^f, v,,^,) = {vf\e'"^*f-^^^\v,), (2.9) 
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where \vi) and \vf) are eigenstates of volume operator in 'Hf™", and cj) is the internal time. As is known, this amplitude equals 
to the physical inner product in the timeless framework [8, 9], i.e., 

A{vf,(f>f; Vi,(j)i) = {vf,(j)f\vi,(f>i)phy, (2.10) 

where \v, (j)) is the eigenstate of both volume operator and scalar operator simultaneously in Hkin, which can be written as 
I w) 1 0) . In this sense, the right hand side of (2. 1 0) can be written as 

/oo 
da{<l>f\{vf\e'"''2\p^\ \v.)\cp,), (2.11) 
-OO 

by inner product (2.6). Sine the constraint operator C has been separated into gravitational part and material part, which live in 
l-L^ln^ and "H]!^^/^ separately, by inserting a set of generalized basis Ip^) of the right hand side of Eq.(2.1 1) becomes 

iaO I 



da{<j,f\e'"^2\p4 / dpM{p^\<j,,){vf\e-'"''\v,) 
= — dp^2bJ e^Tf (^/-^-M dae'"^{vf\e-''"''\v,). (2.12) 

Because we restrict the framework in positive frequency, a unite step function 0{p^) has to be introduced to eUminate negative 
valued Hence we have [7] 

^TTT dp^2\p^\0{p^) e'^^^f-^'U dae'"^{vf\e-''^^\v,). (2.13) 

Thus, the transition amplitude in the deparameterized framework has been formulated by the group averaging. In order to 
compute {vf\e~'''"^\vi), we split the exponential into N identical pieces and insert complete basis, as is done in non-relativistic 
quantum mechanics path integral, 

(i)/|e-'"«|i;,) = {vf\e-''''^\vN-i}{vN-i\e-''''^\vN-2}...{vi\e-''^^\v,}, (2.14) 

DJV-1,...,-U1 

where e = -i-. In the limit oo(e — )• 0), the operator e~""® can be expanded to the second order, and hence we get 

K|e-''^"®|i;„_i) = <5,„„,„„_, - iea{v„\e\v„-i) + ©(e^). (2.15) 

To simplify the calculation, we now employ the Hamiltonian operator 9 = 9o in the simplified LQC of fc = FRW model [11], 
where Qq is a second order difference operator given by 

eo*(w) = v[{v + 2)1'(iJ + 4) - 2v'9{v) + {v- 2)'^{v - 4)]. (2.16) 



This equation leads to 



(Wn|6oK'n-l) = -^-T^tin-i " " ^ (<5t,„ ,i;„_ i +4 " 2Sy,^^y^_^ + (5t,„ ,„„_ i -4 ) ■ (2.17) 



Applying (2.17) to (2.15) and writing Kronecker delta as integral of 6„, which acts as the role of conjugate variable of w„, we 
have 

K|e-'^^"®«|i;„_i) 

=- r db^ e-''"^--— - ^ae^v^^, "-\''-U sin' &„) + 0{e'). 

TT Jo 4 2 

(2.18) 
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Applying (2.14) and (2.18) to (2.13), and taking 'continuum limit', we obtain 

1 nOO nOO p2 

1 r ^ 

lim — Af / dhN-dhi TT 

vi<s -\,...,v\ n— 1 

lim — ^ / db^.-dbi —— / dp02|p0|6'(p0) / da 



X lim > — ^ / dbN-CLb, I I e»6n(fn~t-„-l)/2g-»ae^^„_ r" + ;" ^ 4sin^fe„ 



AT 

X exp i 



e -^(0/ - (pi) - Kyvn ~ Wn-i)/2 + I — ^""-1 ^ 4sm 6„ 



(2.19) 



Supposing that we can integrate out p^, a first and then take the limit, the final result comes up in a path integral formulation as 



|gii?o(0/-0.)|„^) ^ j Vv j Vb cxp^ j d<j){V3^hvsuib-^), (2.20) 



where ii = From Eq.(2.20) we can read out the effective Hamiltonian of the deparameterized system as Tieff = 

—y/iirGhv sin b, which coincides with that from canonical approach to first order 

Let us summarize what we have done to construct the path integral formulation of k=0 LQC. First, the transition amplitude 
is expressed in terms of group averaging as (2.13). Secondly, for deparameterized framework, we cut the gravitational part of 
exponential e~'"® in the group averaging expression into N pieces and then calculate each exponential to the second order Fi- 
nally, we integrate out material variable and the parameter a to get the transition amplitude in terms of exponential integrated 
over all possible paths in gravitational phase space. Note that what we did in the second step is different from that in ordinary 
quantum mechanics which concerns nothing about group averaging. Since our path integral formulation is to be constructed un- 
der deparameterized framework, we regard the material part of e*"*^, i.e., e^°'P<p^^ , as complex function which can be integrated 
out. Hence only the gravitational part e~*"® is split as (2.14). As a matter of fact, one may also split the whole e'""-^ and get 
the path integral formulation under timeless framework [10]. This two kinds of path integral formulation are equivalent to each 
other The two approaches will be compared for spatially curved FRW models in next sections. 



III. PATH INTEGRAL FORMULATION OF fc = 1 LQC 

The spatial manifold S of = 1 FRW model is topologically a 3-sphere S'^. Unlike fc = model, we do not have to 
introduce a fiducial cell V in it. Rather, the volume K> of S with reference to a fiducial metric °qab is finite. We shall set 

-j^/g . (1/3)' 

■— Vo ■ The gravitational connections and the density weighted triads can still be expressed as — cVo °uj\ 
and = pVo \/°q °e°. The fundamental Poisson bracket remains {c, p\ ~ ^i^G^jZ, and we introduce the same new 
conjugate variables as in (2.2). However, the gravitational part of the Hamiltonian constraint of fc = 1 model is different from 
that in fc = model due to non-vanishing spatial curvature. In quantum theory, we have the same kinematical Hilbert space 
^fcin = ® ^fcin"- The gravitational Hamiltonian constraint operator Oi can be written as a combination of diagonal and 

off-diagonal part, whose action reads [17] 



ei*(v) ^ e-'^^-^eoe*^"^ ^(v) + ei*(t;), (3.1) 



where 



/ 4 Ij^l /W, ^ sV 3A ' 

ei = 37rGt;2[- sin^C^) + (1 + 7')(^)2] = Qr(v\ (3.2) 

and 9o is the Hamiltonian operator in fc = model. The definition of K is in accordance with [14]. For simplicity, we 
still employ the simplified treatment and use Eq.(2.16). To calculate the transition amplitude between two basis vectors in the 
deparameterized framework as (2.10), we also need to calculate Eq.(2.15) with 8 replaced by 8i. Note that the diagonal part 
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(vn|6i|wn-i) is easy to deal with by writing delta function as an integration over b. We are going to evaluate the nontrivial 
off-diagonal part 



SttG 



-Vn-l- 



4 " ^ 2 

Since we would like to derive an effective action valid in low energy regime, we assume w ^ 1 and hence get 

J{v + 4) - J{v) = fi{v), f{v - 4) - fiv) = -fi{v). 

Thus, Eq.(3.3) turns out to be 



4 ^" ^ 2 
Consequently, the matrix element of 0i can be calculated as 



1^71-1 / aoe ^ " " 4sm (b - ^(u„_i) — ). 



(3.4) 



(3.5) 



(w,i|9i|w„_i) 
(i6„ e 



4&n(fn-fn-l)/2 



Bi(w„-i) + 37rG u„_i ^ sm (6- /x(u„_i)y) 



(3.6) 



Therefore we have 



-ibn {Vn—Vn-l}/2 



1- iae[ &i{vn-i) + SttG ^ sin^(6 - ^(u„_i)^) 



O(e'). (3.7) 



Applying (3.7) and (2.14) to (2.13), and integrating out and a, finally we arrive at the path integral formulation of /c = 1 
FRW model under deparameterized framework as 



Vv / Vb cxp- / d(?i(ftWei(w) + 37rGu2sin^[5-/l(i;) 



)• 



This enable us to read out an effective action as 



S[bM^ I d4>{mei{v) + 37TGv^sm^[b~fliv)^-^]~^), 



and thus the effective Hamiltonian under deparameterized framework. 



n'^^P = ei(u) + SttG w2 sin' [6 - fi{v)-^ ^. 
Hence the effective equations of this model can be derived by following Poisson brackets 



(3.8) 



(3.9) 



(3.10) 



(3.11a) 



(3.11b) 



where a prime on a letter denotes the derivative with respect to the argument. 

We could also develop path integral formulation of fc = 1 FRW model under timeless framework. For the sake of simplicity, 
group averaging can be re-defined in replacement of (2.6) as 



Gi (w) + 37rG w2 sin^ [b - Ji{v) ^] 



. e;(i;) + 67rGt.sin2[6-/2(tOf ] - iirGv^ sin2[b ~ Ji{vY-^\jJi' {v) 

b - {b, K^ff I . 

ei(w) + 37rG w2 sin2[6-/2(i;)f] 



/Is )vhy 



da 



(3.12) 
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Thus the transition ampUtude can be expressed as 

/•OO _ 

A{vf,(Pf;v^, / o^a(0/|(^;/|e'"^|^',)|(/'^> 



da 



N 



1 E n('^"i(^'"i^""^i^"-i)i'^«-i) 

'Vn — 1 , ■ ■ ■ n— 1 

/ da d(bN-i-.d^i Y[{Me'"'-\'Pn-i){vn\e-'"''''\vn-i), 



: lim / da 

Af— >oo 



(3.13) 



iJiv — n— 1 



where the second step is done by spUtting the whole e'""-^ into pieces and inserting a complete basis. To calculate the material 
part, we insert the basis of jp^^ ) and get 



(3.14) 



Applying (3.7) and (3.14) into (3.13), we obtain 

A{vf,(j)f;vi,(f>,) 



■ lim / da 

W-S-oo 



CJO 



1 



9N^i...a(pi 



{2nh) 



N 



dp^^...dp^^ 



VN-1.---.V1 



N 
ra=l 



+ a I - Oi(u„_i) - SttG U„_ 



dbN---dbi 

Vn + Vn-1 . 2 



2 sin-'(6- /x(u„_i) — ) 



(3.15) 



Taking the 'continuum limit' in which J2 ^ ^ lo '^'''^ finally get the amplitude 

n=l 

A{vf,(t>f;vi,(l)i) 



da I V(f) I Vp^ I Vv I Vb exp - / dr 



hvb 



+ ha['-^~ Qiiv) - SttG sin'^{b - fl{v)^) 



(3.16) 



where a dot over a letter denotes the derivative with respect to t. From this formulation we could read out the effective Hamil- 
tonian constraint of timeless framework as 



tls 



- Qiiv) - SttG sin^ib - fi(v)^) = 



-eff~-f- ^.^.2 

Note that in [17] an effective Hamiltonian constraint under timeless framework was also given by canonical theory as 



(3.17) 



Ceff _ A{V) 



IGttG IGttG 



sm n{c - y ) - sm (— ) + (1 + 7 ) — 



+ (!^)"l5(„)|.0, (3.18) 



where 



A{v) 



27 K Sir I 



6y7jHII--i|-K' + i| 

3 



B{v)= (^) K\v\ |i; + l|i/3_|^,_i|i/3 



(3.19) 



In low energy regime where v ^ 1, we have 



A{v) 



27 K 8tt 



6 73/2 



2v, 



K 



B{v) = — + 0{v- 

V 



(3.20) 
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Then the effective Hamiltonian constraint (3.18) can be simpHfied as 

^'^f V^r^^^i " ^ ^ 

Note that (3.21) is slightly different from (3.17) by a factor ^-^)~^, which will not affect effective equations of motion. 

The appearance of this factor is due to the difference of Eq.(2.3) from Eq.(2.4). The effective canonical equations of the model 
can be derived from (3.21) as 

dv 3^Gz;^sin2[b-/l(z;)^] 

liwy 6 

- = {6, n^s) = . .8...g 3 ' (3.22b) 



Tr = ^"^^ ^ , . ,B-7.-,,| = . ,8...g,.3 • (3.22c) 

Note that we have chosen the plus sign in front of the square-root in Eq.(3.22c), since we focused on 'positive frequency' in 
Eq.(2.8). This means that the inertial time is monotonically increasing with t. It is easy to see that the effective equations 
(3.11) can be derived from equations (3.22) by w = ^^y^ and h ~ ^^y^. Therefore, the effective action in (3.9) predicts the 
same equations of motion for the fc = 1 FRW universe as those of the canonical theory in [17]. The quantum bounce driven by 
the quantum gravity effect would lead to a cyclic universe for this model. Moreover, the comparison between the path integrals 
in the deparameterized framework and timeless framework implies their equivalence. 



e-i = -i?(tO-^ ^"':lw° \v\'l' X II. + 1| - I. - 111. (4.2) 



IV. PATH INTEGRAL FORMULATION OF fc = 1 LQC 

In fc — —1 FRW model, the spatial manifold is noncompact just as that of fc = model, though it is curved. So we should 
still choose a fiducial cell V in the spatial manifold and restrict all integrations on it. The gravitational phase space variables 
A^, E'f and new conjugate variables w, h are defined in the same ways as fc = model. In fc = — 1 LQC [18], by the simplified 
treatment the off diagonal part of the gravitational constraint operator 0_i is exactly the Hamiltonian operator in fc = model, 
while the fc = — 1 model contributes only diagonal part, i.e., 

e_i = eo + e_i, (4.1) 

where 

In the low energy regime (v ^ 1), we have 

e_i(w) = --iT^Gv^-i^v}^^. (4.3) 

In the simplified treatment, we can still apply the 0o in (2.16) to (4.1). Thus, by Eq.(2.17), we can easily write down the matrix 
element of 0_i as 

(v„|e_i|w„_i) 

= i r dh,, e--''"(""-"-^)/^[e,iK,i) + ^„- /" + ''"-^ sin^ /Ic]. (4.4) 

Hence, through the approach outlined in Section 2, the path integral formulation of the deparameterized framework for fc = — 1 
FRW model can be obtained as 

= ^Vv Jvbexp^J d(t){h^ e_i(w) + SttG ^2 sin^ b - (4.5) 
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Thus the effective Hamiltonian can be read out as 



= -hJ 37tGv^ sin' b - j^Vo' p? 



(4.6) 



Though Poisson bracket we obtain the following effective equations 



b = {bXl^} = - 



SttGu^ sin 26 



<^ SnGv^ ^sin' b — ^^V^^ ^' 

QTTGvS\T?b~Q X i^/^TTG-i^V^V^ 

\7rGv^ (^sin'6-72Vo^/i2^ 



(4.7a) 



(4.7b) 



On the other hand, we can also use (4.4) in the path integral formulation of timeless framework as in Section 3 to get the 
transition amplitude of fc = — 1 model as 



= I da I V(t) I Vp^ I Vv I Vb exp - 



dr 



^ + nal^- Q.iiv) - SttG sin^ b 



(4.8) 



The effective Hamiltonian constraint for k 



1 FRW model is then derived as 



SttGw^ (sin^ 6-7^0-' /I^ 



0. 



Also, Eq.(4.9) is in accordance with the Hamiltonian constraint given by the canonical theory in [18], which reads 



(4.9) 



He 



3VH 



■ sin^ (/ic) 



'■i^/\p\yo 



2/3 



- 3ttGv^ (sin^ b - 72Vo^ fi^ 



\P\'/'PM 



9 V / S^TqM j 

"^Kwy 6 



0. 



(4.10) 



It is easy to check that equations (4.7) are equivalent to the effective equations derived by the Hamiltonian (4.9) or (4.10), and 
the classical big bang singularity is also replaced by a quantum bounce in fc = — 1 LQC. 

In conclusion, the path integral formulation for all FRW models of fc = 0, +1, —1 have been obtained. In all these models, one 
can express the transition amplitude in the deparameterized framework in terms of group averaging. The effective Hamiltonian 
in deparameterized framework and the effective Hamiltonian constraint in timeless framework can respectively be derived by the 
path integrals. They lead to same equations of motion and hence coincide with each other The effective Hamiltonian constraints 
given by the canonical theories of the three models are also confirmed by the path integral formulation. Hence the canonical and 
path integral formulations of LQC coincide. 
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that the effective Hamiltonian derived from the path integral in deparameterized framework is equivalent to the 
effective Hamiltonian constraint derived from the path integral in timeless framework, since they lead to same 
equations of motion. Moreover, the effective Hamiltonian constraints of above models derived in canonical 
theory are confirmed by the path integral formulation. 
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I. INTRODUCTION 



A core concept of general relativity is that physical theory should be background independent, which indicates that the gravity 
and matter evolve as a whole entity without anything pre-existed. This idea is being carried out in loop quantum gravity 
(LQG)[l-4], which is a non-perturbative and background independent quantization of general relativity. The construction of 
LQG inspired the research on spin foam models (SFMs) [5], which are proposed as the path-integral formulation for LQG. 
In SFMs, the transition amplitude between physical quantum states is formulated as a sum over histories of all physically 
appropriate states. The heuristic picture is the following. One considers a 4-dimensional spacetime region M bounded by two 
spacelike 3-surfaces, which would correspond to the initial and final states for the path integral. The feature of SFMs is to 
■'nI" , employ the spin network states in LQG to represent the initial and final states. A certain path between the two boundaries is a 
^ quantum 4-geometry interpolated between the two spin networks. The interpolated spin foams can be constructed by considering 
the dual triangulation of M and coloring its surfaces with half integers and edges with suitable intertwiners. In order to obtain 
the physical inner product between the two boundary states, we have to sum over all possible triangulations and colorings [1]. 
\ This is certainly difficult to be achieved. 

On the other hand, the idea and technique of LQG is successfully carried out in the symmetry-reduced models, known as 
Loop Quantum Cosmology (LQC) [6]. Since the infinite degrees of freedom of gravity have been reduced to finite ones by 
homogeneity (and isotropy), LQC provides a simple arena to test the ideas and constructions of the full LQG. Therefore, it 
is also desirable to test the idea and construction of SFMs by LQC models [7]. It was shown in Refs.[8, 9] that the transition 
amplitude in the deparameterized LQC of fc = Friedmann universe equals the physical inner product in the timeless framework, 
. ^ • and concrete evidence has been provided in support of the general paradigm underlying SFMs through the lens of LQC. How 
I to achieve local spinfoam expansion in LQC was also studied in Refs. [10, 11]. Recently, the effective Hamiltonian constraint 
^ . of fc = LQC was derived from the path integral formulation in timeless framework in Ref.[12], and the effective measure 
introduces some conceptual subtleties in arriving at the WKB approximation. The purpose of this paper is to study the path 
integral formulation of LQC mimicking the construction of SFMs and following [7-9, 11, 12]. The new ingredient in this work 
is to formulate the path integrals by multiple group-averaging rather than to divide a single group-averaging. Besides the fc = 
Friedmann-Rorberston- Walker (FRW) model, the fc = ±1 models will also be considered. It turns out that the physical inner 
product between boundary states, defined by group-averaging approach, can be interpreted as sum over all possible quantum 
states. The relation between the path integral formalism in timeless framework and that in deparameterized framework can 
be revealed by the group-averaging. Moreover, we will be able to derive the effective action and thus Hamiltonian of the 
above LQC models by the path integral formalism in timeless framework as well as deparameterized framework. The effective 
Hamiltonian coincide with those obtained by canonical approach. Hence the path integral and canonical formulations of LQC 
confirm each other The effective Hamiltonian constraints of above models in canonical theory are also confirmed by the path 
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integral fonnulation. Starting from the canonical LQC, we will derive the path integral formulation and the effective Hamiltonian 
of fc = 0, +1, —1 FRW models respectively in the following Sections 2, 3, and 4. To simplify the calculations, we will employ 
the simplified LQC treatments [13]. Note that, although the path integrals of fc = FRW LQC have been studied in both 
deparameterized and timeless frameworks, we propose a new formulation by the multiple group-averaging approach. 

II. PATH INTEGRAL FORMULATION OF fc = LQC 

The model that we are considering is the spatially flat (fc — 0) FRW universe filled with a massless scalar field. In the 
kinematical setting, one has to introduce an elementary cell V in the spatial manifold and restricts all integrations to this cell. 
Fix a fiducial flat metric °q^^ and denote by Vo the volume of V in this geometry. The gravitational phase space variables — the 

— fl/3) — f2/3) 

connections and the density weighted triads — can be expressed as = c 14 "w^ and = pVo where 
{°uJa, °ef) are a set of orthonormal co-triads and triads compatible with and adapted to V. p is related to the scale factor a 

via \p\ = Vo^^a^- The fundamental Poisson bracket is given by; {c, p} = 8t:Gj/S, where G is the Newton's constant and 7 
the Barbero-Immirzi parameter The gravitational part of the Hamiltonian constraint reads Cgrav — ~6c^y^|p[/7^. To apply the 
area gap A of full LQG to LQC [14], it is convenient to introduce variable fi satisfying 

fl^\p\ = A=iAV37rj)il (2.1) 

where ip = Gh, and new conjugate variables [13, 15-17]: 

v:=2V3fi^^, b:=flc. (2.2) 

The new canonical pair satisfies {b,v} = j-. On the other hand, the matter phase space consists of canonical variables (j) and 
Pcf, which satisfy = 1. To mimic LQG, the polymer-like representation is employed to quantize the gravity sector The 

kinematical Hilbert space for gravity then reads T-L^]™' = L^(MBohr, d/^if ), where RboHt is the Bhor compactification of the 
real line and d^iH is the Haar measure on it [18]. It turns out that the eigenstates of volume operator v, which are labeled by 
real number v, constitute an orthonomal basis in "H^.^^" as {vi\v2) = ^vi,v2- For the scalar matter sector, one just uses the 
standard Schrodinger representation for its quantization, where the kinematical Hilbert space is W^^^ = i^(R, dip). The total 
kinematical Hilbert Space of the system is a tensor product, T-Lkin = 'Hkln" ® '^ki'n*^ '■^^ above two Hilbert spaces. As a 
totally constrained system, the dynamics of this model is reflected in the Hamiltonian constraint Cgrav + Cmatt = 0. Quantum 
mechanically, physical states are those satisfying quantum constraint equation 

(Cgrav + Cn,att)*(^^, </>) =0. (2.3) 

It is not difficult to write the above equation as [19] 

d^(v, 0) EE (fl - eW(v, (j)) = 0. (2.4) 
On one hand, Eq.(2.4) indicates that we can get physical states by group averaging kinematical states as [10] 

■^f{v,4,) = lirn^ ^ j " da e'"^ /(«,</)), (2.5) 



and thus the physical inner product of two states reads 



f\9 )vhy - ) = lim ^ H da (/|e^"^|5). (2.6) 



As is known, in timeless framework the transition amplitude equals to the physical inner product [8, 9], i.e., 

1 f"" - 

Atis{vf,(t)f;vi,4>i) = {vf,(j)}\ui,(l)^)phy= lim - — / da{vf,4>f\e'°'^\vi,4>i) (2.7) 

tto-^oo Zao J-a„ 

On the other hand, Eq.(2.4) can also be written as 

dl^{v,(j)) + m{v,(t))=Q. (2.8) 
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The similarity between Eq.(2.8) and Klein-Gorden equation suggests that one can regard (f> as internal time, with respect to which 
gravitational field evolves. In this deparameterized framework, we focus on positive frequency solutions, i.e., those satisfying 

~id^^+{v, (j,) - ye*+(v, (j)) = H^+{v, (j,). (2.9) 
The transition amplitude in deparameterized framework is then given by 

Adep{vf,<i>f - tH,^,) = {vf\e'"^'^f-'^'^\v,), (2.10) 

where \vi) and \vf) are eigenstates of volume operator in H^™", and is the internal time. 

Next we are to deal with transition amplitude under timeless work by starting from (2.7). In order to compute 
{vf, (t)f\e^°''~^\vi, (pi), a straightforward way is to split the exponential into N identical pieces and insert complete basis as in 
[12]. However, since a is the group-averaging parameter which goes from — oo to oo, it is unclear whether a could be treated 
as the time variable t in non-relativistic quantum mechanics path integral. We thus consider alternative path integral formulation 
by group-averaging. To this end, we notice the following identity. 



1 r- 

lim / da e 

'2ao J-aa 



iaC\ 



lim -i^ / ° da'-^ / °dae'("'+")^|z;,</)), (2.11) 



which could be generalized as 



1 /""iVo 1 , 

lim / daN...—^ d(3;i e*("i+ -+"")'^|w». (2.12) 



lim — f " dae'°'^\v,(f>) 

a„->oo 2ao J-a„ 

= lim _ .-. 

ajvo,...,0!io^oo 2aNo J-&NC 2^1o 
N 

Let a — ^ <5iT Then this equation indicates that we did not expect each piece of the exponential e'"*^ to be identical. Rather, 

n=l 

they could be independent of each other. In order to trace the power for expansion, we let ci„ = ea„, where e = j^. Then (2.7) 
becomes 

1 1 /■"'" i E 'ia„d 

a,„^oo 2aNo J -UNO ^^'o J -Silo 

= lim / daN--- / dai(u/, (/)/|e "=i I'l'j, (2.13) 

a«„,...,Q,„^oo ZaNo J -UNO ^"■'o J-a,o 

where «„„ = oi„ol^-,n = 1,2, ..,iV. Next, we are going to insert a set of complete basis at each knot. Notice that is 
the eigenstate of both volume operator and scalar operator simultaneously in Hkin, which can be written as \v) |(/)) for short of 
\v) ® 10), and 



Thus, we have 

. " - N 

(«/,</)/|e'"=i'""^|i;„</),) = E j ddpN^i..A^i\{{K\{vnW'''-^K-i)K~i), (2.15) 



vn-\,---V\ n—1 



where Vf ~ VN,(f)f = (f)^, Vi = vo,(j)i = (j)o have been set. Since the constraint operator C has been separated into gravitational 
part and material part, which live in 'H^™'' and "H^""** separately, we could calculate the exponential on each kinematical space 
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separately. For the material part, one gets 



.|e"""^|0„„i) 

- ^1 



(2.16) 



As for the gravitational part, in the limit N oo{e ^ 0), the operator e """^ can be expanded to the first order, and hence we 
get 



(2.17) 



To simplify the calculation, we now employ the Hamiltonian operator = 0o in the simplified LQC of fc = FRW model [13], 
where Oq is a second order difference operator given by 



This equation leads to 



So'^iv) = -^^v[{v + 2)1'(t; + 4) - 2v'i>{v) + {v - 2)1'(w - 4)]. 



(2.18) 



(2.19) 



Applying (2.19) to (2.17) and writing Kronecker delta as integral of 6„, which acts as the role of conjugate variable of w„, we 
have 

K|e-"""§°|z;„-i) 

TT ./n 4 2 



(2.20) 



Applying (2.16) and (2.20) to (2.13), and then taking 'continuum limit', we obtain 



- lim lim 



n 2a„ 



X / d(f>N-l 



i 



1 



N 



daN... / dai 



— OO 

N 



dp4,^...dp^^ e 

n=l 

^ 1 

Y[ exp i 



-4sin^ 6„ 



= lim lim 

AT— >-00 Q«„,...,Qjo— S-OO \ 2a, 



roo 

X / aojv 



duN... I dai 

N 



I dbN--dbi 



N 

X I I expie 

n=l 



l>JV_l,...,tIl 

o2 



Finally, we could write the above equation in path integral formulation as 



Atis{vf,(j)f; Vi,4'i) = c I Va I V4> I Vp^ j Vv j Vh cxp - j dr 



SttG Vn + Vn-l . . 2 

— 77 



4 sin bn 



p,f,(t> ~ —bv + ha{-j^ ~ 3ttGv sin 6) 



(2.21) 



(2.22) 
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where c is certain constant, and a dot over a letter denotes the derivative with respect to r. The effective Hamiltonian constraint 
can be read out from (2.22) as 

Ce f f = ^ - SttGw^ sin^ b, (2.23) 

which coincides with that in Refs.[12, 16]. 

On the other hand, for the deparameterized framework, we may still employ above group-averaging viewpoint. To this end, 
we define a new constraint operator C+ = ^ — H. Then Eq.(2.9) can be rewritten as 

CV*+(w,0)=O. (2.24) 

The transition amplitude for this new constraint reads 



Adep{vf,(j)f;vi,(j)i)^ lim - — / dQ;(v/, 0/|e*"'^+ 0^) = lim - — / da{vf,(f>f\2\p^\9{p^)e'-°''^\vi,(j), 



Qo-i-oo 2ao J -a Qo->tx3 2a, 



If ! 



(2.25) 



where 



\P<t>\\P<t>) = \P<t>\\P<t>) , &{P4>)\P<t>) = \! \ !^ n 
Similar to the timeless case, the integration over single a can be written as multiple integrations, 

N 

= lim / daw---- / dai {vf,4>f\e "=i \vi,4>i). (2.27) 

Since we work in the deparameterized case, it is reasonable to insert in the completeness relation [20, 21] 

l = (2.28) 

V 

where 

= lim — / d/3e'''^+|z;,0). (2.29) 

Note that the Hilbert space in deparameterized framework is unitarily equivalent to the physical Hilbert space with the basis 
(2.29). Then the first piece of the exponential in (2.27) becomes 

= lim / " dai(z;i,,/)i|e""i^|«o>o>, (2.30) 



and the last piece of the exponential reads 

lim / daN^ d/35v(«w,'/'jv|e"""^+e"'«^+|wiv-i,</'jv-i) 

= lim / daw(ww,0w|e"""^+|wAr-i,0w-i). (2.31) 

The remaining pieces of the exponential can also be expressed as 

lim J— f "° dan{Vn,Me'"'-^+\Vn-l,^n-l). (2.32) 



an — >-oo 



2a. 



6 



So (2.27) becomes 

Adep{vf,(j)}; Vi,4>i 



lim 



1 



lim 



da 



N- 



2aio J-ai 



TV 



a«o,...,Q;„— i-oo 2a/Vo 



da 



AT. 



2aio 



E n (^^",</'«|e"""^+|z;„-i, </)„-!). 

^'~)e"""^|i;„_i,0„-i), (2.33) 



^^A^- 1 ,. . . ,^^1 n— 1 



dai ^ (Wn 



|2K„ 



where (2.25) is applied in second step. Now we can split each piece in (2.33) into gravitational and material parts. Calculations 
similar to those in timeless framework lead to 



1 



Adep{vf,(j)f; Vi,(t)i) 
- lim lim 

AT— i-oo aA.„,...,Qj„^oo 2q.No 
N 



X 



n 2b0„|^'(p</,„)cxpie 









pair, 




dajy. 






2a]„ 


' —C^Io 




K - <? 




b„ V 


h 


e 




2 



dai 



1 

27rh 



N 




SttG 



1>JV_1,...,1>1 

V„ + Vn-1 



dbN-.-dbi 



4 sin 6i, 



(2.34) 



We can integrate out a„ and and arrive at 



^ 1 



: lim lim 

\ri=l 

N 

X expie{(j)f - 0;) 

n=l 




-Vn-1 



2 eicl^f - 



Vv I Vb cxp — / d<j){V SttGHv sin b 



hvb. 
'~2~' 



(2.35) 



where v = The effective Hamiltonian in deparameterized framework can be read out from (2.35) as -ffe// = VSttGIIv sin b, 
which is just the deparameterized Hamiltonian of (2.23). 

Let us summarize what we have done to construct the path integral formulation of k=0 LQC. First, the transition amplitude 
is expressed in terms of multiple group averaging. Secondly, we insert in complete basis and calculate each exponential to 
the second order Thirdly, by taking 'continuum limit', the path integral formulation under timeless framework is derived. 
This technique could also be applied to the path integral calculation of other constrained systems. Finally, by reexpressing the 
deparameterized equation in the form of constraint, the path integral of deparameterized framework is also formulated in this 
approach. Hence we have got path integral formulation of timeless framework as well as deparameterized framework. The 
relation of these two kinds of path integral formulation is revealed by (2.25). They are in accordance with each other in the sense 
of effective theory. 

III. PATH INTEGRAL FORMULATION OF fc = 1 LQC 

The spatial manifold M of fc = 1 FRW model is topologically a 3-sphere S'^. Unlike fc = model, we do not have to 
introduce a fiducial cell in it. Rather, the volume Vo of M with reference to a fiducial metric °qab is finite. We shall set 
'■= Vo ■ The gravitational connections and the density weighted triads can still be expressed as = cVo 

and i?f = pVo y/°q °e°. The fundamental Poisson bracket remains {c, p} ~ &tiG^ /i, and we introduce the same new 
conjugate variables as in (2.2). However, the gravitational part of the Hamiltonian constraint of fc = 1 model is different from 
that in fc = model due to non-vanishing spatial curvature. In quantum theory, we have the same Hilbert space T-Lkin ~ 
'^kin' ® ^fein"- The gravitational Hamiltonian constraint operator 8i can be written as a combination of diagonal and off- 
diagonal part, whose action reads [22] 



(3.1) 
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where 



/ 4 Ij^l Jy^)^ ^ 3 V 3A ' 

ei = STTG^^h sin^(^) + (1 + 7^)(^)^] ^ e,(.), 



(3.2) 



and 00 = is the Hamiltonian operator in fc = model. The definition of K is in accordance with [17]. For simplicity, 
we still employ the simplified treatment and use Eq.(2.18). To calculate the transition amplitude between two basis vectors 
in the timeless framework as (2.7), we also need to calculate Eq.(2.17) with 8 replaced by 9i. Note that the diagonal part 
(v„|8i|w„_i) is easy to deal with by writing delta function as an integration over b. We are going to evaluate the nontrivial 
off-diagonal part. 



SttG 



-Vn-l- 



Vn + Vn-1 



4 ""^ 2 
We assume w > 4 and do Taylor expansion as 

w2/3 /2 4 1 2, 1_4 



f{v + 4)-f{v) = 



4a \6v2\66v 3! 3 3 6 v 



fiv - 4) - f{v) 



Ik 



2 4 12 14 12 1 4 4 
6v 2! 3 6 V 3! 3 3 6 v 



-g{v) + h{v), 



(3.4) 



where h{v) denotes sum of the odd terms and g{v) that of the evens. Thus, Eq.(3.3) turns out to be 



SttG Vn + Vn-l 1 

Vn-l 



db„ e 



-ihn {Vn—Vn-l)/2 



g-^foM..-i)2cos(26„ - - 2 



4 2 

If we denote A„ ~ e~*^°''(^"-i)2 cos(2&„ — (7(w„_i)^o) ~ 2, the matrix element of Oi can be calculated as 



(3.5) 



(w„|8i|u„_i 
1 



c(6„ e " 8i(u„_i)H ^ w„_i A, 



(3.6) 



Therefore we have 



TT 



\-iae\ 8i(w„_i) H — v^-i A, 



+ 0(e2). 



(3.7) 



Replacing (2.20) with (3.7) and following the procedure as in fc = model, finally we arrive at the path integral formulation of 
k = \ FRW model under timeless framework as 



k=l 



ci Va V(j) Vp^ Vv Vb cxp 



^ Jo 



dr 



hvb . I ~ 3ttG 2 . / , N 
- Oi(v) — V A{v,b) 



ha 



(3.8) 



However, since A(z;, b) is a complex function, (3.8) would imply a complex effective action or Hamiltonian. To avoid this 
problem, we notice that in semiclassical regime, the higher power terms in (3.4) decrease dramatically as v increases, and hence 
only the first term, denoted by gi, dominates. For example, when v = 100, ^^^'^^^^(■^lYiv) ''^^ < 1%. As a matter of fact, in the 
simplified LQC treatment a similar approximation has already been used. In this sense, we could reserve only gi term and get 



A{v,b) = 2 [cos{2b - gi{v)£o) - 1] ^ Ash-i^ {b - fl{v)£o/2). 



(3.9) 
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Consequently, the transition amplitude turns out to be 

Atls{vf,(f>f;V^,(j)i)''=^ 



ci Va V(t) Vp^ Vv Vb exp 



dr 



hvb 



ha 



h 



61(1;) -SttG siii^ib- fliv)^) 



(3.10) 



where a dot over a letter denotes the derivative with respect to t. From this formulation we could read out the effective Hamil- 
tonian constraint of timeless framework as 



c: 



lis 



n 



Qi(v) - SttG siii2(& - ^^(v) — ). 



Note that in [22] an effective Hamiltonian constraint under timeless framework was also given by canonical theory as 



"He 



C. 



A{v) 



16nG 16ttG 



sm ^(c - y ) - sm ( — ) + (1 + 7 



(3.11) 



(3.12) 



where 



A{v) 
B{v) 



27K fSi^ £p . ... . . 



-I KM 



\v + l\ 



1/3 



\v-l 



1/3 



(3.13) 



In large v regime, we have 



27K /Stt ip „ 



6 73/2 



Then the effective Hamiltonian constraint (3.12) can be simplified as 

Hel 



(3.14) 



(3.15) 



Note that (3.15) is slightly different from (3. 11) by a factor which will not affect effective equations of motion. 

The appearance of this factor is due to the difference of Eq.(2.3) from Eq.(2.4). 

Taking into account of Heff = 0, the effective canonical equations can be derived from (3.15) as 



dv 
d^ 

db 
d7 



dr 



{v,Heff} - 

{b,neff}-- 



SttGv^ sin2[6-^(w)f ] 



KK 



^\ 6 I 



e'l [v) + GttG V sin^ [b - fl{v)^] - SttGw^ sin 2[b - p,{v)^]fl'{v) 



"ffP'V 6 ' ' 



61 (w) + SttG i;2sin^[6 - Tl{v)^] 



^ 1 6 ' Kh^ \ 6 r 



(3.16a) 
(3.16b) 

(3.16c) 



where a prime on a letter denotes the derivative with respect to the argument. Note that in order to compare the effective 
equations with those in deparameterized framework, we focused on 'positive frequency' as Eq.(2.9) and thus chose the plus sign 
in front of the square-root in Eq.(3. 16c). This means that the internal time is monotonically increasing with r. 

Following the same procedure in /s = model, we also obtain the path integral formulation of k = 1 model under deparame- 
terized framework as 



/ Vv / Vb exp 



\fc=i 



d<p{\lQi{v) + 37rG z;2 sin2[6 - fl{v)^]h - ^) 



(3.17) 
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This enable us to read out an effective action as 



2rj -/ \^0i^ hvb.^ 



S[b,v]= I d(b{\lQiiv) + 3TTGv^sin'[b-fl{v)^]h-—), (3.18) 



and thus the effective Hamiltonian under deparameterized framework, 



H^^P = -^Qi{v) + 37TGv^sm^[b~fliv)^-^]h. (3.19) 
Hence the effective equations of this model can be derived by following Poisson brackets 

— ^{v,n^/}^ (3.20a) 

^/ei(w) + 3^G«2sin2[6-/2(i;)|] 



_ ud'^PX - e; jv) + 6nG V sin^ [b - fljvyf] - inGv' sin 2 [6 - /2(z;) |]/2- (z;) 



^ei{v) + SnG sm^[b - ] 



It is easy to see that the effective equations (3.20) can be derived from equations (3.16) by = d^^/dr ^"'^ ^ ~ d^/dr ■ 
Therefore, the effective action in (3.18) predicts the same equations of motion for the k = 1 FRW universe as those in [22]. The 
quantum bounce driven by the quantum gravity effect would lead to a cyclic universe for this model. Moreover, the comparison 
between the path integrals in the deparameterized framework and timeless framework implies their equivalence in the sense of 
effective theory. 

IV. PATH INTEGRAL FORMULATION OF fc = 1 LQC 

In fc = — 1 FRW model, the spatial manifold is noncompact just as that of fc = model. So we should still choose a fiducial 
cell V in the spatial manifold and restrict all integrations on it. The gravitational phase space variables (AJ^, Ef) and reduced 
new conjugate variables (v, b) are defined in the same ways as fc = model. In fc = — 1 LQC [23], by the simplified treatment 
the off diagonal part of the gravitational constraint operator 0_i is exactly the Hamiltonian operator in fc = model, while the 
fc = — 1 model contributes only diagonal part, i.e. 

e_i = eo + e_i, (4.1) 

where 

12x1/3 

In large v regime, we have 

e_i(w) = -SttGuVVq^M^- (4.3) 

In the simpHfied treatment, we can still apply the 60 in (2.18) to (4.1). Thus, by Eq.(2.19), we can easily write down the matrix 
element of 0_i as 

{Vn\Q-l\Vn-l) 

= - r db^ e-'^"(""-'"-^)/^[e_iK_i) + 3^Gz;„_ /""^^''"~^ sinVc]. (4.4) 
Hence, through the approach outlined in Section 2, we can also obtain the path integral formulation of timeless framework as 

1 " / 2 \ " 



e_i = -i3(t,)-i^^^^A_|„|i/3||„ + i| _ |„ _ i||. (4.2) 



=c_i j Va \ V(j) f Vp^ I Vv I Vb cxp - / dr 







2 ' '""^ \ h 



^ + na{^~ inGv^ ( sin^ b - fi^ 
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where a dot over a letter denotes the derivative with respect to r. The effective Hamiltonian constraint for fc = — 1 FRW model 
is then derived as 

2 

%%■ - ^ - ZttGv^ (sin^ b - i''v}fi^) . (4.6) 
Also, Eq.(4.6) is in accordance with the Hamiltonian constraint given in [23], which reads 

a2/3 



^ ■ 2,- , ^ ^ I ,3/2 



^-3^Gt;2 (sin^ 6-72^^5/22) 



(4.7) 



On the other hand, we also use (4.4) to obtain the path integral formulation of A: = — 1 model under deparameterized frame- 
work as 



Adep{vf,(j)f\Vi,(j)i) 



k=-l 



=c'_i j Vv j Vh exp ^ j d(j){^J 3ttGv^ (sin^ b ~ -fWj fl^^ h - ^). (4.8) 
The effective Hamiltonian in the deparameterized framework for k = —1 FRW model can be derived as 

W^f = -^3ttGv^ (sin2 b ~ jW^^p^^h. (4.9) 
Though Poisson bracket we obtain the following effective equations 

r -lydepi 37rGu2sin26 m ^ 

V = {v,n^/} = ^ =, (4.10a) 

JsnGv^ (sin^ b - jW^^ fi^^ 
b = {b,n^J} = 6.G.sin^5-6x3V3.G,X^,i ^^^^^^ 
SttGv^ (sin^ b - j^o^ fL^ 

It is easy to check that equations (4.10) are equivalent to the effective equations derived by the Hamiltonian (4.6) or (4.7), and 
the classical big bang singularity is also replaced by a quantum bounce in fc = —1 LQC. 

In conclusion, a multiple group-averaging approach is proposed to derive the path integral formulation for all FRW models of 
fc = 0, +1, —1. In all these models, one can express the transition amplitude in the deparameterized framework in terms of group 
averaging. The relation between the path integral formulation in timeless framework and that in deparameterized framework is 
then clarified. The effective Hamiltonian in deparameterized framework and the effective Hamiltonian constraint in timeless 
framework can respectively be derived by the path integrals. They lead to same equations of motion and hence coincide with 
each other The effective Hamiltonian constraints given by the canonical theories of the three models are also confirmed by the 
path integral formulation. Hence the canonical and path integral formulations of LQC coincide. 
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We study the path integral formulation of Friedmann universe filled with a massless scalar field in loop 
quantum cosmology. All the isotropic models of fc = 0, +1, —1 are considered. To construct the path integrals 
in the timeless framework, a multiple group- averaging approach is proposed. Meanwhile, since the transition 
amplitude in the deparameterized framework can be expressed in terms of group-averaging, the path integrals 
can be formulated for both deparameterized and timeless frameworks. Their relation is clarified. It turns out 
that the effective Hamiltonian derived from the path integral in deparameterized framework is equivalent to the 
effective Hamiltonian constraint derived from the path integral in timeless framework, since they lead to same 
equations of motion. Moreover, the effective Hamiltonian constraints of above models derived in canonical 
theory are confirmed by the path integral formulation. 
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I. INTRODUCTION 



A core concept of general relativity is that physical theory should be background independent, which indicates that the gravity 
and matter evolve as a whole entity without anything pre-existed. This idea is being carried out in loop quantum gravity 
(LQG)[l-4], which is a non-perturbative and background independent quantization of general relativity. The construction of 
LQG inspired the research on spin foam models (SFMs) [5], which are proposed as the path-integral formulation for LQG. 
In SFMs, the transition amplitude between physical quantum states is formulated as a sum over histories of all physically 
appropriate states. The heuristic picture is the following. One considers a 4-dimensional spacetime region M bounded by two 
spacelike 3-surfaces, which would correspond to the initial and final states for the path integral. The feature of SFMs is to 
■'nI" , employ the spin network states in LQG to represent the initial and final states. A certain path between the two boundaries is a 
^ quantum 4-geometry interpolated between the two spin networks. The interpolated spin foams can be constructed by considering 
the dual triangulation of M and coloring its surfaces with half integers and edges with suitable intertwiners. In order to obtain 
the physical inner product between the two boundary states, we have to sum over all possible triangulations and colorings [1]. 
\ This is certainly difficult to be achieved. 

On the other hand, the idea and technique of LQG is successfully carried out in the symmetry-reduced models, known as 
Loop Quantum Cosmology (LQC) [6]. Since the infinite degrees of freedom of gravity have been reduced to finite ones by 
homogeneity (and isotropy), LQC provides a simple arena to test the ideas and constructions of the full LQG. Therefore, it 
is also desirable to test the idea and construction of SFMs by LQC models [7]. It was shown in Refs.[8, 9] that the transition 
amplitude in the deparameterized LQC of fc = Friedmann universe equals the physical inner product in the timeless framework, 
. ^ • and concrete evidence has been provided in support of the general paradigm underlying SFMs through the lens of LQC. How 
I to achieve local spinfoam expansion in LQC was also studied in Refs. [10, 11]. Recently, the effective Hamiltonian constraint 
^ . of fc = LQC was derived from the path integral formulation in timeless framework in Ref.[12], and the effective measure 
introduces some conceptual subtleties in arriving at the WKB approximation. The purpose of this paper is to study the path 
integral formulation of LQC mimicking the construction of SFMs and following [7-9, 11, 12]. The new ingredient in this work 
is to formulate the path integrals by multiple group-averaging rather than to divide a single group-averaging. Besides the fc = 
Friedmann-Rorberston- Walker (FRW) model, the fc = ±1 models will also be considered. It turns out that the physical inner 
product between boundary states, defined by group-averaging approach, can be interpreted as sum over all possible quantum 
states. The relation between the path integral formalism in timeless framework and that in deparameterized framework can 
be revealed by the group-averaging. Moreover, we will be able to derive the effective action and thus Hamiltonian of the 
above LQC models by the path integral formalism in timeless framework as well as deparameterized framework. The effective 
Hamiltonian coincide with those obtained by canonical approach. Hence the path integral and canonical formulations of LQC 
confirm each other The effective Hamiltonian constraints of above models in canonical theory are also confirmed by the path 
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integral fonnulation. Starting from the canonical LQC, we will derive the path integral formulation and the effective Hamiltonian 
of fc = 0, +1, —1 FRW models respectively in the following Sections 2, 3, and 4. To simplify the calculations, we will employ 
the simplified LQC treatments [13]. Note that, although the path integrals of fc = FRW LQC have been studied in both 
deparameterized and timeless frameworks, we propose a new formulation by the multiple group-averaging approach. 

II. PATH INTEGRAL FORMULATION OF fc = LQC 

The model that we are considering is the spatially flat (fc — 0) FRW universe filled with a massless scalar field. In the 
kinematical setting, one has to introduce an elementary cell V in the spatial manifold and restricts all integrations to this cell. 
Fix a fiducial flat metric °q^^ and denote by Vo the volume of V in this geometry. The gravitational phase space variables — the 

— fl/3) — f2/3) 

connections and the density weighted triads — can be expressed as = c 14 "w^ and = pVo where 
{°uJa, °ef) are a set of orthonormal co-triads and triads compatible with and adapted to V. p is related to the scale factor a 

via \p\ = Vo^^a^- The fundamental Poisson bracket is given by; {c, p} = 8t:Gj/S, where G is the Newton's constant and 7 
the Barbero-Immirzi parameter The gravitational part of the Hamiltonian constraint reads Cgrav — ~6c^y^|p[/7^. To apply the 
area gap A of full LQG to LQC [14], it is convenient to introduce variable fi satisfying 

fl^\p\ = A=iAV37rj)il (2.1) 

where ip = Gh, and new conjugate variables [13, 15-17]: 

v:=2V3fi^^, b:=flc. (2.2) 

The new canonical pair satisfies {b,v} = j-. On the other hand, the matter phase space consists of canonical variables (j) and 
Pcf, which satisfy = 1. To mimic LQG, the polymer-like representation is employed to quantize the gravity sector The 

kinematical Hilbert space for gravity then reads T-L^]™' = L^(MBohr, d/^if ), where RboHt is the Bhor compactification of the 
real line and d^iH is the Haar measure on it [18]. It turns out that the eigenstates of volume operator v, which are labeled by 
real number v, constitute an orthonomal basis in "H^.^^" as {vi\v2) = ^vi,v2- For the scalar matter sector, one just uses the 
standard Schrodinger representation for its quantization, where the kinematical Hilbert space is W^^^ = i^(R, dip). The total 
kinematical Hilbert Space of the system is a tensor product, T-Lkin = 'Hkln" ® '^ki'n*^ '■^^ above two Hilbert spaces. As a 
totally constrained system, the dynamics of this model is reflected in the Hamiltonian constraint Cgrav + Cmatt = 0. Quantum 
mechanically, physical states are those satisfying quantum constraint equation 

(Cgrav + Cn,att)*(^^, </>) =0. (2.3) 

It is not difficult to write the above equation as [19] 

d^(v, 0) EE (fl - eW(v, (j)) = 0. (2.4) 
On one hand, Eq.(2.4) indicates that we can get physical states by group averaging kinematical states as [10] 

■^f{v,4,) = lirn^ ^ j " da e'"^ /(«,</)), (2.5) 



and thus the physical inner product of two states reads 



f\9 )vhy - ) = lim ^ H da (/|e^"^|5). (2.6) 



As is known, in timeless framework the transition amplitude equals to the physical inner product [8, 9], i.e., 

1 f"" - 

Atis{vf,(t)f;vi,4>i) = {vf,(j)}\ui,(l)^)phy= lim - — / da{vf,4>f\e'°'^\vi,4>i) (2.7) 

tto-^oo Zao J-a„ 

On the other hand, Eq.(2.4) can also be written as 

dl^{v,(j)) + m{v,(t))=Q. (2.8) 
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The similarity between Eq.(2.8) and Klein-Gorden equation suggests that one can regard (f> as internal time, with respect to which 
gravitational field evolves. In this deparameterized framework, we focus on positive frequency solutions, i.e., those satisfying 

~id^^+{v, (j,) - ye*+(v, (j)) = H^+{v, (j,). (2.9) 
The transition amplitude in deparameterized framework is then given by 

Adep{vf,<i>f - tH,^,) = {vf\e'"^'^f-'^'^\v,), (2.10) 

where \vi) and \vf) are eigenstates of volume operator in H^™", and is the internal time. 

Next we are to deal with transition amplitude under timeless work by starting from (2.7). In order to compute 
{vf, (t)f\e^°''~^\vi, (pi), a straightforward way is to split the exponential into N identical pieces and insert complete basis as in 
[12]. However, since a is the group-averaging parameter which goes from — oo to oo, it is unclear whether a could be treated 
as the time variable t in non-relativistic quantum mechanics path integral. We thus consider alternative path integral formulation 
by group-averaging. To this end, we notice the following identity. 



1 r- 

lim / da e 

'2ao J-aa 



iaC\ 



lim -i^ / ° da'-^ / °dae'("'+")^|z;,</)), (2.11) 



which could be generalized as 



1 /""iVo 1 , 

lim / daN...—^ d(3;i e*("i+ -+"")'^|w». (2.12) 



lim — f " dae'°'^\v,(f>) 

a„->oo 2ao J-a„ 

= lim _ .-. 

ajvo,...,0!io^oo 2aNo J-&NC 2^1o 
N 

Let a — ^ <5iT Then this equation indicates that we did not expect each piece of the exponential e'"*^ to be identical. Rather, 

n=l 

they could be independent of each other. In order to trace the power for expansion, we let ci„ = ea„, where e = j^. Then (2.7) 
becomes 

1 1 /■"'" i E 'ia„d 

a,„^oo 2aNo J -UNO ^^'o J -Silo 

= lim / daN--- / dai(u/, (/)/|e "=i I'l'j, (2.13) 

a«„,...,Q,„^oo ZaNo J -UNO ^"■'o J-a,o 

where «„„ = oi„ol^-,n = 1,2, ..,iV. Next, we are going to insert a set of complete basis at each knot. Notice that is 
the eigenstate of both volume operator and scalar operator simultaneously in Hkin, which can be written as \v) |(/)) for short of 
\v) ® 10), and 



Thus, we have 

. " - N 

(«/,</)/|e'"=i'""^|i;„</),) = E j ddpN^i..A^i\{{K\{vnW'''-^K-i)K~i), (2.15) 



vn-\,---V\ n—1 



where Vf ~ VN,(f)f = (f)^, Vi = vo,(j)i = (j)o have been set. Since the constraint operator C has been separated into gravitational 
part and material part, which live in 'H^™'' and "H^""** separately, we could calculate the exponential on each kinematical space 
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separately. For the material part, one gets 



.|e"""^|0„„i) 

- ^1 



(2.16) 



As for the gravitational part, in the limit N oo{e ^ 0), the operator e """^ can be expanded to the first order, and hence we 
get 



(2.17) 



To simplify the calculation, we now employ the Hamiltonian operator = 0o in the simplified LQC of fc = FRW model [13], 
where Oq is a second order difference operator given by 



This equation leads to 



So'^iv) = -^^v[{v + 2)1'(t; + 4) - 2v'i>{v) + {v - 2)1'(w - 4)]. 



(2.18) 



(2.19) 



Applying (2.19) to (2.17) and writing Kronecker delta as integral of 6„, which acts as the role of conjugate variable of w„, we 
have 

K|e-"""§°|z;„-i) 

TT ./n 4 2 



(2.20) 



Applying (2.16) and (2.20) to (2.13), and then taking 'continuum limit', we obtain 



- lim lim 



n 2a„ 



X / d(f>N-l 



i 



1 



N 



daN... / dai 



— OO 

N 



dp4,^...dp^^ e 

n=l 

^ 1 

Y[ exp i 



-4sin^ 6„ 



= lim lim 

AT— >-00 Q«„,...,Qjo— S-OO \ 2a, 



roo 

X / aojv 



duN... I dai 

N 



I dbN--dbi 



N 

X I I expie 

n=l 



l>JV_l,...,tIl 

o2 



Finally, we could write the above equation in path integral formulation as 



Atis{vf,(j)f; Vi,4'i) = c I Va I V4> I Vp^ j Vv j Vh cxp - j dr 



SttG Vn + Vn-l . . 2 

— 77 



4 sin bn 



p,f,(t> ~ —bv + ha{-j^ ~ 3ttGv sin 6) 



(2.21) 



(2.22) 
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where c is certain constant, and a dot over a letter denotes the derivative with respect to r. The effective Hamiltonian constraint 
can be read out from (2.22) as 

Ce f f = ^ - SttGw^ sin^ b, (2.23) 

which coincides with that in Refs.[12, 16]. 

On the other hand, for the deparameterized framework, we may still employ above group-averaging viewpoint. To this end, 
we define a new constraint operator C+ = ^ — H. Then Eq.(2.9) can be rewritten as 

CV*+(w,0)=O. (2.24) 

The transition amplitude for this new constraint reads 



Adep{vf,(j)f;vi,(j)i)^ lim - — / dQ;(v/, 0/|e*"'^+ 0^) = lim - — / da{vf,(f>f\2\p^\9{p^)e'-°''^\vi,(j), 



Qo-i-oo 2ao J -a Qo->tx3 2a, 



If ! 



(2.25) 



where 



\P<t>\\P<t>) = \P<t>\\P<t>) , &{P4>)\P<t>) = \! \ !^ n 
Similar to the timeless case, the integration over single a can be written as multiple integrations, 

N 

= lim / daw---- / dai {vf,4>f\e "=i \vi,4>i). (2.27) 

Since we work in the deparameterized case, it is reasonable to insert in the completeness relation [20, 21] 

l = (2.28) 

V 

where 

= lim — / d/3e'''^+|z;,0). (2.29) 

Note that the Hilbert space in deparameterized framework is unitarily equivalent to the physical Hilbert space with the basis 
(2.29). Then the first piece of the exponential in (2.27) becomes 

= lim / " dai(z;i,,/)i|e""i^|«o>o>, (2.30) 



and the last piece of the exponential reads 

lim / daN^ d/35v(«w,'/'jv|e"""^+e"'«^+|wiv-i,</'jv-i) 

= lim / daw(ww,0w|e"""^+|wAr-i,0w-i). (2.31) 

The remaining pieces of the exponential can also be expressed as 

lim J— f "° dan{Vn,Me'"'-^+\Vn-l,^n-l). (2.32) 



an — >-oo 



2a. 
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So (2.27) becomes 

Adep{vf,(j)}; Vi,4>i 



lim 



1 



lim 



da 



N- 



2aio J-ai 



TV 



a«o,...,Q;„— i-oo 2a/Vo 



da 



AT. 



2aio 



E n (^^",</'«|e"""^+|z;„-i, </)„-!). 

^'~)e"""^|i;„_i,0„-i), (2.33) 



^^A^- 1 ,. . . ,^^1 n— 1 



dai ^ (Wn 



|2K„ 



where (2.25) is applied in second step. Now we can split each piece in (2.33) into gravitational and material parts. Calculations 
similar to those in timeless framework lead to 



1 



Adep{vf,(j)f; Vi,(t)i) 
- lim lim 

AT— i-oo aA.„,...,Qj„^oo 2q.No 
N 



X 



n 2b0„|^'(p</,„)cxpie 









pair, 




dajy. 






2a]„ 


' —C^Io 




K - <? 




b„ V 


h 


e 




2 



dai 



1 

27rh 



N 




SttG 



1>JV_1,...,1>1 

V„ + Vn-1 



dbN-.-dbi 



4 sin 6i, 



(2.34) 



We can integrate out a„ and and arrive at 



^ 1 



: lim lim 

\ri=l 

N 

X expie{(j)f - 0;) 

n=l 




-Vn-1 



2 eicl^f - 



Vv I Vb cxp — / d<j){V SttGHv sin b 



hvb. 
'~2~' 



(2.35) 



where v = The effective Hamiltonian in deparameterized framework can be read out from (2.35) as -ffe// = VSttGIIv sin b, 
which is just the deparameterized Hamiltonian of (2.23). 

Let us summarize what we have done to construct the path integral formulation of k=0 LQC. First, the transition amplitude 
is expressed in terms of multiple group averaging. Secondly, we insert in complete basis and calculate each exponential to 
the second order Thirdly, by taking 'continuum limit', the path integral formulation under timeless framework is derived. 
This technique could also be applied to the path integral calculation of other constrained systems. Finally, by reexpressing the 
deparameterized equation in the form of constraint, the path integral of deparameterized framework is also formulated in this 
approach. Hence we have got path integral formulation of timeless framework as well as deparameterized framework. The 
relation of these two kinds of path integral formulation is revealed by (2.25). They are in accordance with each other in the sense 
of effective theory. 

III. PATH INTEGRAL FORMULATION OF fc = 1 LQC 

The spatial manifold M of fc = 1 FRW model is topologically a 3-sphere S'^. Unlike fc = model, we do not have to 
introduce a fiducial cell in it. Rather, the volume Vo of M with reference to a fiducial metric °qab is finite. We shall set 
'■= Vo ■ The gravitational connections and the density weighted triads can still be expressed as = cVo 

and i?f = pVo y/°q °e°. The fundamental Poisson bracket remains {c, p} ~ &tiG^ /i, and we introduce the same new 
conjugate variables as in (2.2). However, the gravitational part of the Hamiltonian constraint of fc = 1 model is different from 
that in fc = model due to non-vanishing spatial curvature. In quantum theory, we have the same Hilbert space T-Lkin ~ 
'^kin' ® ^fein"- The gravitational Hamiltonian constraint operator 8i can be written as a combination of diagonal and off- 
diagonal part, whose action reads [22] 



(3.1) 
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where 



/ 4 Ij^l Jy^)^ ^ 3 V 3A ' 

ei = STTG^^h sin^(^) + (1 + 7^)(^)^] ^ e,(.), 



(3.2) 



and 00 = is the Hamiltonian operator in fc = model. The definition of K is in accordance with [17]. For simplicity, 
we still employ the simplified treatment and use Eq.(2.18). To calculate the transition amplitude between two basis vectors 
in the timeless framework as (2.7), we also need to calculate Eq.(2.17) with 8 replaced by 9i. Note that the diagonal part 
(v„|8i|w„_i) is easy to deal with by writing delta function as an integration over b. We are going to evaluate the nontrivial 
off-diagonal part. 



SttG 



-Vn-l- 



Vn + Vn-1 



4 ""^ 2 
We assume w > 4 and do Taylor expansion as 

w2/3 /2 4 1 2, 1_4 



f{v + 4)-f{v) = 



4a \6v2\66v 3! 3 3 6 v 



fiv - 4) - f{v) 



Ik 



2 4 12 14 12 1 4 4 
6v 2! 3 6 V 3! 3 3 6 v 



-g{v) + h{v), 



(3.4) 



where h{v) denotes sum of the odd terms and g{v) that of the evens. Thus, Eq.(3.3) turns out to be 



SttG Vn + Vn-l 1 

Vn-l 



db„ e 



-ihn {Vn—Vn-l)/2 



g-^foM..-i)2cos(26„ - - 2 



4 2 

If we denote A„ ~ e~*^°''(^"-i)2 cos(2&„ — (7(w„_i)^o) ~ 2, the matrix element of Oi can be calculated as 



(3.5) 



(w„|8i|u„_i 
1 



c(6„ e " 8i(u„_i)H ^ w„_i A, 



(3.6) 



Therefore we have 



TT 



\-iae\ 8i(w„_i) H — v^-i A, 



+ 0(e2). 



(3.7) 



Replacing (2.20) with (3.7) and following the procedure as in fc = model, finally we arrive at the path integral formulation of 
k = \ FRW model under timeless framework as 



k=l 



ci Va V(j) Vp^ Vv Vb cxp 



^ Jo 



dr 



hvb . I ~ 3ttG 2 . / , N 
- Oi(v) — V A{v,b) 



ha 



(3.8) 



However, since A(z;, b) is a complex function, (3.8) would imply a complex effective action or Hamiltonian. To avoid this 
problem, we notice that in semiclassical regime, the higher power terms in (3.4) decrease dramatically as v increases, and hence 
only the first term, denoted by gi, dominates. For example, when v = 100, ^^^'^^^^(■^lYiv) ''^^ < 1%. As a matter of fact, in the 
simplified LQC treatment a similar approximation has already been used. In this sense, we could reserve only gi term and get 



A{v,b) = 2 [cos{2b - gi{v)£o) - 1] ^ Ash-i^ {b - fl{v)£o/2). 



(3.9) 
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Consequently, the transition amplitude turns out to be 

Atls{vf,(f>f;V^,(j)i)''=^ 



ci Va V(t) Vp^ Vv Vb exp 



dr 



hvb 



ha 



h 



61(1;) -SttG siii^ib- fliv)^) 



(3.10) 



where a dot over a letter denotes the derivative with respect to t. From this formulation we could read out the effective Hamil- 
tonian constraint of timeless framework as 



c: 



lis 



n 



Qi(v) - SttG siii2(& - ^^(v) — ). 



Note that in [22] an effective Hamiltonian constraint under timeless framework was also given by canonical theory as 



"He 



C. 



A{v) 



16nG 16ttG 



sm ^(c - y ) - sm ( — ) + (1 + 7 



(3.11) 



(3.12) 



where 



A{v) 
B{v) 



27K fSi^ £p . ... . . 



-I KM 



\v + l\ 



1/3 



\v-l 



1/3 



(3.13) 



In large v regime, we have 



27K /Stt ip „ 



6 73/2 



Then the effective Hamiltonian constraint (3.12) can be simplified as 

Hel 



(3.14) 



(3.15) 



Note that (3.15) is slightly different from (3. 11) by a factor which will not affect effective equations of motion. 

The appearance of this factor is due to the difference of Eq.(2.3) from Eq.(2.4). 

Taking into account of Heff = 0, the effective canonical equations can be derived from (3.15) as 



dv 
d^ 

db 
d7 



dr 



{v,Heff} - 

{b,neff}-- 



SttGv^ sin2[6-^(w)f ] 



KK 



^\ 6 I 



e'l [v) + GttG V sin^ [b - fl{v)^] - SttGw^ sin 2[b - p,{v)^]fl'{v) 



"ffP'V 6 ' ' 



61 (w) + SttG i;2sin^[6 - Tl{v)^] 



^ 1 6 ' Kh^ \ 6 r 



(3.16a) 
(3.16b) 

(3.16c) 



where a prime on a letter denotes the derivative with respect to the argument. Note that in order to compare the effective 
equations with those in deparameterized framework, we focused on 'positive frequency' as Eq.(2.9) and thus chose the plus sign 
in front of the square-root in Eq.(3. 16c). This means that the internal time is monotonically increasing with r. 

Following the same procedure in /s = model, we also obtain the path integral formulation of k = 1 model under deparame- 
terized framework as 



/ Vv / Vb exp 



\fc=i 



d<p{\lQi{v) + 37rG z;2 sin2[6 - fl{v)^]h - ^) 



(3.17) 
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This enable us to read out an effective action as 



2rj -/ \^0i^ hvb.^ 



S[b,v]= I d(b{\lQiiv) + 3TTGv^sin'[b-fl{v)^]h-—), (3.18) 



and thus the effective Hamiltonian under deparameterized framework, 



H^^P = -^Qi{v) + 37TGv^sm^[b~fliv)^-^]h. (3.19) 
Hence the effective equations of this model can be derived by following Poisson brackets 

— ^{v,n^/}^ (3.20a) 

^/ei(w) + 3^G«2sin2[6-/2(i;)|] 



_ ud'^PX - e; jv) + 6nG V sin^ [b - fljvyf] - inGv' sin 2 [6 - /2(z;) |]/2- (z;) 



^ei{v) + SnG sm^[b - ] 



It is easy to see that the effective equations (3.20) can be derived from equations (3.16) by = d^^/dr ^"'^ ^ ~ d^/dr ■ 
Therefore, the effective action in (3.18) predicts the same equations of motion for the k = 1 FRW universe as those in [22]. The 
quantum bounce driven by the quantum gravity effect would lead to a cyclic universe for this model. Moreover, the comparison 
between the path integrals in the deparameterized framework and timeless framework implies their equivalence in the sense of 
effective theory. 

IV. PATH INTEGRAL FORMULATION OF fc = 1 LQC 

In fc = — 1 FRW model, the spatial manifold is noncompact just as that of fc = model. So we should still choose a fiducial 
cell V in the spatial manifold and restrict all integrations on it. The gravitational phase space variables (AJ^, Ef) and reduced 
new conjugate variables (v, b) are defined in the same ways as fc = model. In fc = — 1 LQC [23], by the simplified treatment 
the off diagonal part of the gravitational constraint operator 0_i is exactly the Hamiltonian operator in fc = model, while the 
fc = — 1 model contributes only diagonal part, i.e. 

e_i = eo + e_i, (4.1) 

where 

12x1/3 

In large v regime, we have 

e_i(w) = -SttGuVVq^M^- (4.3) 

In the simpHfied treatment, we can still apply the 60 in (2.18) to (4.1). Thus, by Eq.(2.19), we can easily write down the matrix 
element of 0_i as 

{Vn\Q-l\Vn-l) 

= - r db^ e-'^"(""-'"-^)/^[e_iK_i) + 3^Gz;„_ /""^^''"~^ sinVc]. (4.4) 
Hence, through the approach outlined in Section 2, we can also obtain the path integral formulation of timeless framework as 

1 " / 2 \ " 



e_i = -i3(t,)-i^^^^A_|„|i/3||„ + i| _ |„ _ i||. (4.2) 



=c_i j Va \ V(j) f Vp^ I Vv I Vb cxp - / dr 







2 ' '""^ \ h 



^ + na{^~ inGv^ ( sin^ b - fi^ 
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where a dot over a letter denotes the derivative with respect to r. The effective Hamiltonian constraint for fc = — 1 FRW model 
is then derived as 

2 

%%■ - ^ - ZttGv^ (sin^ b - i''v}fi^) . (4.6) 
Also, Eq.(4.6) is in accordance with the Hamiltonian constraint given in [23], which reads 

a2/3 



^ ■ 2,- , ^ ^ I ,3/2 



^-3^Gt;2 (sin^ 6-72^^5/22) 



(4.7) 



On the other hand, we also use (4.4) to obtain the path integral formulation of A: = — 1 model under deparameterized frame- 
work as 



Adep{vf,(j)f\Vi,(j)i) 



k=-l 



=c'_i j Vv j Vh exp ^ j d(j){^J 3ttGv^ (sin^ b ~ -fWj fl^^ h - ^). (4.8) 
The effective Hamiltonian in the deparameterized framework for k = —1 FRW model can be derived as 

W^f = -^3ttGv^ (sin2 b ~ jW^^p^^h. (4.9) 
Though Poisson bracket we obtain the following effective equations 

r -lydepi 37rGu2sin26 m ^ 

V = {v,n^/} = ^ =, (4.10a) 

JsnGv^ (sin^ b - jW^^ fi^^ 
b = {b,n^J} = 6.G.sin^5-6x3V3.G,X^,i ^^^^^^ 
SttGv^ (sin^ b - j^o^ fL^ 

It is easy to check that equations (4.10) are equivalent to the effective equations derived by the Hamiltonian (4.6) or (4.7), and 
the classical big bang singularity is also replaced by a quantum bounce in fc = —1 LQC. 

In conclusion, a multiple group-averaging approach is proposed to derive the path integral formulation for all FRW models of 
fc = 0, +1, —1. In all these models, one can express the transition amplitude in the deparameterized framework in terms of group 
averaging. The relation between the path integral formulation in timeless framework and that in deparameterized framework is 
then clarified. The effective Hamiltonian in deparameterized framework and the effective Hamiltonian constraint in timeless 
framework can respectively be derived by the path integrals. They lead to same equations of motion and hence coincide with 
each other The effective Hamiltonian constraints given by the canonical theories of the three models are also confirmed by the 
path integral formulation. Hence the canonical and path integral formulations of LQC coincide. 
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